Dark-Energy Dynamics Required to Solve the Cosmic Coincidence 
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Dynamic dark energy (DDE) models are often designed to solve the cosmic coincidence (why, just 
now, is the dark energy density pde, the same order of magnitude as the matter density p m ?) by 
guaranteeing pde ~ pm for significant fractions of the age of the universe. However, such behaviour is 
neither sufficient nor necessary to solve the coincidence problem. Cosmological processes constrain 
the epochs during which observers can exist. Therefore, what must be shown is that a significant 
fraction of observers see pde ~ pm- Precisely when, and for how long, must a DDE model have 
pde ~ pm in order to solve the coincidence? We explore the coincidence problem in dynamic dark 
energy models using the temporal distribution of terrestrial-planet-bound observers. We find that 
any realistic DDE model which can be parameterized as w = wq + w a (l — a) over a few e-folds, 
has pde ~ Pm for a significant fraction of observers in the universe. This demotivates DDE models 
specifically designed to solve the coincidence using long or repeated periods of pde ~ pm- 



I. INTRODUCTION 

In 1998, using supernovae la as standard candles, Ricss 
et al. [l| and Perlmutter et al. @ revealed a recent and 
continuing epoch of cosmic acceleration - strong evidence 
that Einstein's cosmological constant A, or something 
else with comparable negative pressure pde ~ — Pde, cur- 
rently dominates the energy density of the universe ||. 
A is usually interpreted as the energy of zero-point quan- 
tum fluctuations in the vacuum [j, Q with a constant 
equation of state w = p del Pde = — 1. This necessary ad- 
ditional energy component, construed as A or otherwise, 
has become generically known as "dark energy" (DE). 

A plethora of observations have been used to constrain 
the free parameters of the new standard cosmological 
model, ACDM, in which A does play the role of the dark 
energy. Hinshaw et al. [6j] find that the universe is ex- 
panding at a rate of Hq = 71 ± 4 km/s/Mpc; that it 
is spatially flat and therefore critically dense (Cl to to = 
Pi^s- = ^Mptoto = 1-01 ± 0.01); and that the total den- 

PcritO ^ 

sity is comprised of contributions from vacuum energy 
(Cl A0 = 0.74 ± 0.02), cold dark matter (CDM; CIcdmo = 
0.22 ± 0.02), baryonic matter (Cl b0 = 0.044 ± 0.003) and 
radiation (Cl r0 = 4.5 ± 0.2 x 10~ 5 ). Henceforth we will 
assume that the universe is flat {Citato = 1) as predicted 
by inflation and supported by observations. 

Two problems have been influential in moulding ideas 
about dark energy, specifically in driving interest in al- 
ternatives to ACDM. The first of these problems is con- 
cerned with the smallness of the dark energy density 
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[H Hi • Despite representing more than 70% of the total 
energy of the universe, the current dark energy density is 
~ 120 orders of magnitude smaller than energy scales at 
the end of inflation (or ~ 80 orders of magnitude smaller 
than energy scales at the end of inflation if this occurred 
at the GUT rather than Planck scale) 0]. Dark energy 
candidates are thus challenged to explain why the ob- 
served DE density is so small. The standard idea, that 
the dark energy is the energy of zero-point quantum fluc- 
tuations in the true vacuum, seems to offer no solution 
to this problem. 

The second cosmological constant problem 0, [l(| [H| is 
concerned with the near coincidence between the current 
cosmological matter density {p m o ~ 0.26 x p cr ito) and 
the dark energy density (pdeo ~ 0.74 x p CT ito)- In the 
standard ACDM model, the cosmological window during 
which these components have comparable density is short 
(just 1.5 e-folds of the cosmological scalefactor a) since 
matter density dilutes as p m oc a~ 3 while vacuum density 
Pde is constant |12j. Thus, even if one explains why the 
DE density is much less than the Planck density (the 
smallness problem) one must explain why we happen to 
live during the time when pd e ~ p m - 

The likelihood of this coincidence depends on the range 
of times during which we suppose we might have lived. 
In works addressing the smallness problem, Weinberg 
EH considered a multiverse consisting of a large 
number of big bangs, each with a different value of pde- 
There he asked, suppose that we could have arisen in 
any one of these universes; What value of pde should we 
expect our universe to have? While Weinberg supposed 
we could have arisen in another universe, we are simply 
supposing that we could have arisen in another time in 
our universe. We ask, what time t b s , and correspond- 
ing densities Pde(tobs) and p m {tobs) should we expect to 
observe? Weinberg's key realization was that not every 
universe was equally probable: those with smaller pd e 
contain more Milky- Way-like galaxies and are therefore 



time t [s] 



more hospitable @, [H)]. Subsequently, he, and other au- 
thors used the relative number of Milky- Way-like galax- 
ies to estimate the distribution of observers as a func- 
tion of pde, and determined that our value of pde was 
indeed likely [H [H, 

Our value of pd e could have 
been found to be unlikely and this would have ruled out 
the type of multiverse being considered. Here we apply 
analogous reasoning to the cosmic coincidence problem. 
Our observerhood could not have happened at any time 
with equal probability [l2l ]. By estimating the temporal 
distribution of observers we can determine whether the 
observation of pde ~ Pm was likely. If we find pde ~ Pm 
to be unlikely while considering a particular DE model, 
that will enable us to rule out that DE model. 

In a previous paper [12] , we tested ACDM in this way 
and found that pde ~ Pm is expected. In the present 
paper we apply this test to dynamic dark energy models 
to see what dynamics is required to solve the coincidence 
problem when the temporal distribution of observers is 
being considered. 

The smallness of the dark energy density has been an- 
thropically explained in multiverse models with the ar- 
gument that in universes with much larger DE compo- 
nents, DE driven acceleration starts earlier and precludes 
the formation of galaxies and large scale structure. Such 
universes are probably devoid of observers [H, [l6|, [l7| • A 
solution to the coincidence problem in this scenario was 
outlined by Garriga et al. [18| who showed that if pd e is 
low enough to allow galaxies to form, then observers in 
those galaxies will observe r ~ 1. 

To quantify the time-dependent proximity of p m and 
Pde, we define a proximity parameter, 



Pde Pn 
Pm ' Pd 



(1) 



which ranges from r ~ 0, when many orders of magni- 
tude separate the two densities, to r = 1, when the two 
densities are equal. The presently observed value of this 
parameter is tq — ~ 0.35. In terms of r, the coinci- 
dence problem is as follows. If we naively presume that 
the time of our observation t i>s has been drawn from 
a distribution of times Pt{t) spanning many decades of 
cosmic scalefactor, we find that the expected proximity 
parameter is r « <C 0.35. In the top panel of Fig. Q] we 
use a naive distribution for t f, s that is constant in log(a) 
to illustrate how observing r as large as ro ~ 0.35 seems 
unexpected. 

In Lineweaver and Egan we showed how the appar- 
ent severity of the coincidence problem strongly depends 
upon the distribution Pt(t) from which t t, s is hypoth- 
esized to have been drawn. Naive priors for t bs, such 
as the one illustrated in the top panel of Fig. [U lead to 
naive conclusions. Following the reasoning of Weinberg 
0) 0) 03 we interpret Pt(t) as the temporal distribu- 
tion of observers. The temporal and spatial distribution 
of observers has been estimated using large (lO n M ) 
galaxies [l3|, HE HI, OH and terrestrial planets [l2[ as 
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FIG. 1: (Top) The history of the energy density of the uni- 
verse according to standard ACDM. The dotted line shows 
the energy density in radiation (photons, neutrinos and other 
relativistic modes). The radiation density dilutes as a -4 as 
the universe expands. The dashed line shows the density in 
ordinary non-relativistic matter, which dilutes as a -3 . The 
thick solid line shows the energy of the vacuum (the cosmo- 
logical constant) which has remained constant since the end 
of inflation. The thin solid peaked curve shows the proximity 
r of the matter density to the vacuum energy density (see 
Eq. [T}. The proximity r is only ~ 1 for a brief period in the 
log(a) history of the cosmos. Whether or not there is a coin- 
cidence problem depends on the distribution Pt(t) for t b s - If 
one naively assumes that we could have observed any epoch 
with equal probability (the light grey shade) then we should 
not expect to observe r as large as we do. If, however, Pt(t) is 
based on an estimate of the temporal distribution of observers 
(the dark grey shade) then ro ~ 0.35 is not surprising, and the 
coincidence problem is solved under ACDM [12§. (Bottom) 
The dark energy density history is modified in DDE models. 
Observational constraints on the dark energy density history 
are represented by the light grey shade (details in Section lrTT)) . 



tracers. The top panel of Fig. [T] shows the temporal dis- 
tribution of observers Pt(t) from Lineweaver and Egan 
0. 

A possible extension of the concordance cosmologi- 
cal model that may explain the observed smallness of 
Pde is the generalization of dark energy candidates to 
include dynamic dark energy (DDE) models such as 
quintessence, phantom dark energy, k-essence and Chap- 
lygin gas. In these models the dark energy is treated as 
a new matter field which is approximately homogenous, 
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and evolves as the universe expands. DDE evolution of- 
fers a mechanism for the decay of p deify from the ex- 
pected Planck scales (10 93 g/cm 3 ) in the early universe 
(1CP 44 s) to the small value we observe today (10~ 30 
g/cm 3 ). The light grey shade in the bottom panel of 
Fig- [H represents contemporary observational constraints 
on the DDE density history. Many DDE models are 
designed to solve the coincidence problem by having 
p deify ~ Pmify for a large fraction of the history /future 
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of the universe [IE 

Pde ~ Pm for extended or repeated periods the hope is to 
ensure that r ~ 1 is expected. 



Our main goal in this paper is to take into account the 
temporal distribution of observers to determine when, 
and for how long, a DDE model must have pd e ~ Pm 
in order to solve the coincidence problem? Specifically, 
we extend the work of Lineweaver and Egan [12J to find 
out for which cosmologies (in addition to ACDM) the 
coincidence problem is solved when the temporal distri- 
bution of observers is considered. In doing this we answer 
the question, Does a dark energy model fitting contem- 
porary constraints on the density pd e and the equation 
of state parameters, necessarily solve the cosmic coinci- 
dence? Both positive and negative answers have inter- 
esting consequences. An answer in the affirmative will 
simplify considerations that go into DDE modeling: any 
DDE model in agreement with current cosmological con- 
straints has pde ~ Pm for a significant fraction of ob- 
servers. An answer in the negative would yield a new 
opportunity to constrain the DE equation of state pa- 
rameters more strongly than contemporary cosmological 
surveys. 



II. DYNAMIC DARK ENERGY MODELS IN 
THE FACE OF THE COSMIC COINCIDENCE 

Though it is beyond the scope of this article to pro- 
vide a complete review of DDE (see Copeland et al. 
[46| . Szydlowski et al. [131), here we give a few representa- 
tive examples in order to set the context and motivation 
of our work. Fig. [2] illustrates density histories typical 
of tracker quintessence, tracking oscillating energy, in- 
teracting quintessence, phantom dark energy, k-essence, 
and Chaplygin gas. They are discussed in turn below. 



A. Quintessence 

In quintessence models the dark energy is interpreted 
as a homogenous scalar field with Lagrangian density 
£(cj), X) = ±</> 2 - V{(t>) [3 HI M HH M, M 0- The 
evolution of the quintessence field and of the cosmos de- 
pends on the postulated potential V((f>) of the field and on 
any postulated interactions. In general, quintessence has 

a time- varying equation of state w = — — ^Jt—tMt- 

J b H Pde <t> 2 /2 + V{<)>) 

Since the kinetic term <j> 2 /2 cannot be negative, the equa- 
tion of state is restricted to values w > — 1. Moreover, 
if the potential V(<p) is non-negative then w is also re- 
stricted to values w < +1. 

If the quintessence field only interacts gravitationally 
then energy density evolves as = — 3(w + 1)^ and 
the restrictions — 1 < w < +1 mean pd e decays (but 
never faster than c 
increases) . 



or remains constant (but never 



Tracker Quintessence 



A different coincidence problem arises when the time 
of observation is conditioned on and the parameters of a 
model are allowed to slide. The tuning of parameters and 
the necessity to include ad-hoc physics are large problems 
for many current dark energy models. This paper does 
not address such issues, and the interested reader is re- 
ferred to Hebecker and Wetterich [43| , Bludman [44| and 
Linder .45]. In the coincidence problem addressed here 
we let the time of observation vary to see if r(t bs) > 0.35 
is unlikely according to the model. 

In Section|TT]we present several examples of DDE mod- 
els used to solve the coincidence problem. An overview of 
observational constraints on DDE is given in Section [Till 
In Section IIVI we estimate the temporal distribution of 
observers. Our main analysis is presented in Section [V] 
Our main result - that the coincidence problem is solved 
for all DDE models fitting observational constraints - is 
illustrated in Fig. [7] Finally, in Sect ion fVTI we end with a 
discussion of our results, their implications and potential 
caveats. 



Particular choices for V(tf>) lead to interesting attrac- 
tor solutions which can be exploited to make pde scale 
("track") sub-dominantly with p r + p m . 

The DE can be forced to transit to a A- like (w w —1) 
state at any time by fine-tuning V(<p). In the A-like 
state pde overtakes p m and dominates the recent and 
future energy density of the universe. We illustrate 
tracker quintessence in Fi g. [2] using a power law potential 
V{4>) — M(j)~ a (panel b) |49l. l5ll. |53]| and an exponential 
potential V(<p) = Mexp(l/0) (panel c) [20[. 

The tracker paths are attractor solutions of the equa- 
tions governing the evolution of the field. If the tracker 
quintessence field is initially endowed with a density off 
the tracker path (e.g. an equipartition of the energy avail- 
able at reheating) its density quickly approaches and 
joins the tracker solution. 

2. Oscillating Dark Energy 

Dodelson et al. [2(| explored a quintessence 
potential with oscillatory perturbations V{4>) = 
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FIG. 2: The energy density history of the universe according to ACDM (panel a), and seven DDE models selected from the 
literature (see text for references). In each panel the radiation and matter densities are the dotted and dashed lines respectively. 
The DE density is given by the thick black line. The proximity parameter r is given by the thin black line at the base of each 
panel. Of the DDE models shown here, tracker quintessence and k-essence (panels b, c and g) have r ~ 1 for a small fraction of 
the life of the universe (whether the abscissa is t, log(i), a, log(a), or any other of a large number of measures). On the other 
hand, tracking oscillating energy, interacting quintessence, phantom DE and Chaplygin gas (panels d, e, f and h) exhibit r ~ 1 
for a large fraction of the life of the universe. For the phantom DE example (panel f) this is true in t, but not in a or log(a). In 
phantom models the future universe grows super-exponentially to a = oo (a "big-rip") shortly after matter-DE equality. Thus 
the universe spends a large fraction of time with r ~ 1, however this is is not seen in log(a)-space. For each of the models in 
this figure, numerical values for free parameters were chosen to crudely fit observational constraints and are given in Appendix 

El 



Mexp(— X4>) [1 + Asm(i/<f>)]. They refer to models of 
this type as tracking oscillating energy. Without the 
perturbations (setting A = 0) this potential causes exact 
tracker behaviour: the quintessence energy decays as 
Pr + Pm and never dominates. With the perturbations 



the quintessence energy density oscillates about p r + p m 
as it decays (Fig. - The quintessence energy dom- 
inates on multiple occasions and its equation of state 
varies continuously between positive and negative values. 
One of the main motivations for tracking oscillating 
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energy is to solve the coincidence problem by ensuring 
that pde ~ Pm or pd e ~ p r at many times in the past or 
future. 

It has yet to be seen how such a potential might arise 
from particle physics. Phcnomenologically similar cos- 
mologies have been discussed in Ahmed et al. (26|, Feng 
et al. [13, Yang and Wang 



large (12%). Note that for the "naive Pf,(t)" prior shown 
in Fig. [TJ the solution of Caldwell et al. [6l[ and Scherrer 
[33| fails because r > 0.1 is brief in log(a)-space. It fails 
in this way for many other choices of Pt (t) including, for 
example, distributions constant in a or log(t). 

C. K-Essence 



3. Interacting Quintessence 

Non-gravitational interactions between the 
quintessence field and matter fields might allow en- 
ergy to transfer between these components. Such 
interactions are not forbidden by any known symmetry 
[56| . The primary motivation for the exploration 
of interacting dark energy models is to solve the 
coincidence problem. In these models the present mat- 
ter/dark energy densit y p roximity r may be constant 

di m, pi pa nn m ii m, m m [Mm or s iowi y 

varying [29l. |35||. 

We plot a density history of the interacting 
quintessence model of Amendola [Hj] in Fig. This 
model is characterized by a DE potential V(4>) — 
Aexp[Bc/)\ and DE-matter interaction term Q = —Cp m <p, 
specifying the rate at which energy is transferred to the 
matter fields. The free parameters were tuned such 
that radiation domination ends at a = 1CP 5 and that 
rt^oo — 0.35. 



In k-essence the DE is modeled as a scalar field with 
non-canonical kinetic energy [g^, HH, HH, (65|. Non- 
canonical kinetic terms can arise in the effective ac- 
tion of fields in string and supergravity theories. Fig. 
[5^ shows a density history typical of k-essence models. 
This particular model is from Armcndariz-Picon et al. 
[64| and Steinhardt [ill ]. During radiation domination 
the k-essence field tracks radiation sub-dominantly (with 
Wde — w r — 1/3) as do some of the other models in 
Fig. O However, no stable tracker solution exists for 
Wde — w m {— 0). Thus after radiation-matter equality, 
the field is unable to continue tracking the dominant 
component, and is driven to another attractor solution 
(which is generically A-like with Wd e ~ — 1). The onset 
of DE domination was recent in k-essence models because 
matter-radiation equality prompts the transition to a A- 
like state. K-essence thereby avoids fine-tuning in any 
particular numerical parameters, but the Lagrangian has 
been constructed ad-hoc. 



B. Phantom Dark Energy 

The analyses of Riess et al. [H|] and Wood-Vasey et al. 
[5^ 1 have mildly (~ la) favored a dark energy equation of 
state Wde < — I- These values are unattainable by stan- 
dard quintessence models but can occur in phantom dark 
energy models [60], in which kinetic energies are negative. 
The energy density in the phantom field increases with 
scalefactor, typically leading to a future "big rip" sin- 
gularity where the scalefactor becomes infinite in finite 
time. Fig. [5f shows the density history of a simple phan- 
tom model with a constant equation of state w = —1.25. 
The big rip (a = 00 at t = 57.5 Gyrs) is not seen in 
log(a)-space. 



Caldwell et al. [6l| and Scherrer [33j have explored 



how phantom models may solve the coincidence problem: 
since the big rip is triggered by the onset of DE domi- 
nation, such cosmologies spend a significant fraction of 
their total time with r large. For the phantom model with 
w = -1.25 (Fig. [5f) Scherrer [H finds r > 0.1 for 12% 
of the total lifetime of the universe. Whether this solves 
the coincidence or not depends upon the prior probability 
distribution Pt(t) for the time of observation. Caldwell 
et al. [6l| and Scherrer [33[ implicitly assume that the 
temporal distribution of observers is constant in time (i.e. 
Pt(t) = constant). For this prior the coincidence prob- 
lem is solved because the chance of observing r > 0.1 is 



D. Chaplygin Gas 

A special fluid known as Chaplygin gas motivated by 
braneworld cosmology may be able to play the role of 
dark matter and the dark energy [66l . |67T ] . Generalized 
Chaplygin gas has the equation of state pd e = —^Pde 
which behaves like pressureless dark matter at early times 
(wde ~ when pd e is large), and like vacuum energy at 
late times {wde ~ — 1 when pd e is small). In Fig. [2h we 
show an example with a = 1. 



E. Summary of DDE Models 

Two broad classes of DDE models emerge from our 
comparison: 

1. In ACDM, tracker quintessence and k-essence mod- 
els, the dark energy density is vastly different from 
the matter density for most of the lifetime of the 
universe (panels a, b, c, g of Fig. [2]). The coinci- 
dence problem can only be solved if the probability 
distribution P t (t) for the time of observation is nar- 
row, and overlaps significantly with anr~l peak. 
If Pt(t) is wide, e.g. constant over the life of the uni- 
verse in t or log(t), then observing r ~ 1 would be 
unlikely in these models and the coincidence prob- 
lem is not resolved. 



6 



2. Tracking oscillating energy, interacting 
quintessence, phantom models and Chaplygin 
gas models (panels d, e, f, h of Fig. [2]) employ 
mechanisms to ensure that r ~ 1 for large fractions 
of the life of the universe. In these models the 
coincidence problem may be solved for a wider 
range of Pt(t) including, depending on the DE 
model, distributions that are constant over the 
whole life of the universe in t, log(t), a or log(a). 

The importance of an estimate of the distribution P t (t) is 
highlighted: such an estimate will either rule out models 
of the first category because they do not solve the coin- 
cidence problem, or demotivate models of the second be- 
cause their mechanisms are unnecessary to solve the co- 
incidence problem. This analysis does not addressed the 
problems associated with fine-tuning, initial conditions 
or ad hoc mechanisms of many DDE models [HI, [3 [H| . 

We leave this line of enquiry temporarily to discuss 
contemporary observational constraints on the dark en- 
ergy density history, because we wish to test what DE 
dynamics are required to solve the coincidence, beyond 
those which models must exhibit to satisfy standard cos- 
mological observations. 



III. CURRENT OBSERVATIONAL 
CONSTRAINTS ON DYNAMIC DARK ENERGY 

A. Supernovae la 

Observationally, possible dark energy dynamics is ex- 
plored almost solely using measurements of the cosmic 
expansion history. Recent cosmic expansion is directly 
probed by using type la supernova (SNIa) as standard 
candles [l|, 0] • Each observed SNIa provides an indepen- 
dent measurement of the luminosity distance di to the 
redshift of the supernova z$n- The luminosity distance 
to z$n is given by 



where 



E(z) 



di(z S N) = (1 + zsn)- 



H{z) 
H 



N dz 
H J Z=0 Ejzj 



(2) 



(3) 



firop. + zf + n ro0 (i + zf + n de0 



Pde(z) 



PdeO 



and thus depends on Ho, f2 m o, and the evolution of the 
dark energy Pde{z) / Pdeo- The radiation term, irrelevant 
at low redshifts, can be dropped from Equation[3j £l de o is 
a dependent parameter due to flatness (Q de o — 1 — Omo)- 
Contemporary datasets include ~ 200 supernovae at red- 
shifts zsn < 2.16 (a > 0.316) [H, [69| and provide an 
effective continuum of constraints on the expansion his- 
tory over that range [zO, [7l| . The redshift range probed 
by SNIa is indicated in both panels of Fig. [3] 



B. Cosmic Microwave Background 

The first peak in the cosmic microwave background 
(CMB) temperature power spectrum corresponds to den- 
sity fluctuations on the scale of the sound horizon at the 
time of recombination. Subsequent peaks correspond to 
higher- frequency harmonics. The locations of these peaks 
in Z-space depend on the comoving scale of the sound 
horizon at recombination, and the angular distance to re- 
combination. This is summarized by the so-called CMB 
shift parameter R [zl, [73[ which is related to the cosmol- 
ogy by 



R = y/n 



mO 



Zr.e dz 
= W) 



(4) 



where 



1089 7J1 is the redshift of recombina- 



tion. The 3-year WMAP data gives a shift parameter 
R = 1.71 ± 0.03 (El, El- Since the dependence of Equa- 
tion [4] on Hq and fl m o differs from that of Equation [2j 
measurements of the CMB shift parameter can be used 
to break degeneracies between Ho, £l m o an d DE evolu- 
tion in the analysis of SNIa. In the top panel of Fig. [3] we 
represent the CMB observations using a bar from z = 
to z rec . 



Baryonic Acoustic Oscillations and Large Scale 
Structure 



As they imprinted acoustic peaks in the CMB, the 
baryonic oscillations at recombination were expected to 
leave signature wiggles - baryonic acoustic oscillations 
(BAO) - in the power spectrum of galaxies [75[. These 
were detected with significant confidence in the SDSS lu- 
minous red galaxy power spectrum 76]. The expected 
BAO scale depends on the scale of the sound horizon at 
recombination, and on transverse and radial scales at the 
mean redshift zbao, of galaxies in the survey. Eisenstein 
et al. [76|] measured the quantity 



A(z BAO ) 



mO 



E(ZBAO)'- 



1 



zbao 



dz 

W) 



(•5) 



to have a value A(z BAO = 0.35) = 0.469 ± 0.017, thus 
constraining the matter density and the dark energy evo- 
lution parameters in a configuration which is complemen- 
tary to the CMB shift parameter and the SNIa luminos- 
ity distance relation. Ongoing BAO projects have been 
designed specifically to produce stronger constraints on 
the dark energy equation of state parameter w. For ex- 
ample, WIGGLEZ [7tJ will use a sample of high-redshift 
galaxies to measure the BAO scale at zbao ~ 0.75. As 
well as reducing the effects of non-linear clustering, this 
redshift is at a larger angular distance, making the ob- 
served scale more sensitive to w. Constraints from the 
BAO scale depend on the evolution of the universe from 
z rec to zbao to set the physical scale of the oscillations. 
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FIG. 3: The energy densities of radiation p r , matter p m and the cosmological constant pA axe shown as a function of scalefactor, 
by the dotted, dashed, and solid lines respectively. Cosmological probes of dark energy include SNIa, CMB, BAO, the LSS 
linear growth factor and constraints from BBN (see text). Each of these probes is sensitive to the effects of dark energy over 
different redshift intervals, as indicated. The light grey band envelopes wo-w a -parameterized DDE models allowed at < 2a by 
Davis et al. [68l ] (the contour in wo — w a space is shown explicitly in F ig. \7§ . The dark grey band envelopes wo-parameterized 
DDE models (w a — assumed) allowed at < 2a by Wood-Vasey et al. [5S2|. The constraint is w = —1.09 ± 0.16 at 2a. 



They also depend on the evolution of the universe from 
zbao to z = 0, since the observed angular extent of the 
oscillations depends on this evolution. The bar repre- 
senting BAO scale observations in the top panel of Fig. 
[3] indicates both these regimes. 

The amplitude of the B AOs - the amplitude of the large 
scale structure (LSS) power spectrum - is determined 
by the amplitude of the power spectrum at recombina- 
tion, and how much those fluctuations have grown (the 



transfer function) between z rec and zbao- By comparing 
the recombination power spectrum (from CMB) with the 
galaxy power spectrum, the LSS linear growth factor can 
be measured and used to constrain the expansion history 
of the universe (independently of the BAO scale) over 
this redshift range. In practice, biases hinder precise nor- 
malization of the galaxy power spectrum, weakening this 
technique. The range over which this technique probes 
the DE is indicated in Fig. [3] 
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D. Ages 

Cosmological parameters from SNla, CMB, LSS, BAO 
and other probes allow us to calculate the current age of 
the universe to be 13.8±0.1 @ assuming ACDM. Uncer- 
tainties on the age calculated in this way grow dramati- 
cally if we drop the assumption that the DE is vacuum 
energy (w = —1). 

An independent lower limit on the current age of the 
universe is found by estimating the ages of the oldest 
known globular clusters [78| . These observations rule out 
models which predict the universe to be younger than 
12.7 ± 0.7 Gyrs (2a confidence): 

> 12.7 ±0.7 Gyrs. 

Other objects can also be used to set this age limit 
Lineweaver [79|, but generally less successfully due to 
uncertainties in dating techniques. 

Assuming ACDM, an age of 12.7 Gyrs corresponds to a 
redshift of z ~ 5.5. Contemporary age measurements are 
sensitive to the dark energy content from z w 5.5 to z = 
0. In the top panel of Fig.[3]we show this redshift interval. 
The evolution and energy content of the universe before 
12.7 Gyrs ago is not probed by these age constraints. 



zation is w(a) — wq +w a (J — a) [82l. |83|. which does not 
diverge at high redshift. 

Linder and Huterer have argued that the ex- 
tension of this approach to second order, e.g. w(a) = 
Wq + w a (l — a) + w aa (l — a) 2 , is not motivated by current 
DDE models. Moreover, they have shown that next gen- 
eration observations are unlikely to be able to distinguish 
the quadratic from a linear expansion of w. Riess et al. 
[69l | have illustrated this recently using new SNla. 

An alternative technique for exploring the history of 
dark energy is to constrain w(z) or pde(z) in indepen- 
dent redshift bins. This technique makes fewer assump- 
tions about the specific shape of w(z). In the absence 
of any strongly motivated parameterization of w(z) this 
bin-wise method serves as a good reminder of how lit- 
tle we actually know from observation. Using luminos- 
ity distance measurements from SNla, DE evolution has 
been constrained in this way in Az ~ 0.5 bins out to 
redshift zsv ~ 2 [o9l l85l |86{ ■ In the future, BAO mea- 
surements at various redshifts may contribute to these 
constraints, however zbao will probably never be larger 
than zsn- Moreover, because the recombination redshift 
z rec sa 1089 is fixed, only the cumulative effect (from 
z = z rec to z = 0) of the DE can be measured with the 
CMB and LSS linear growth factor. With only this single 
data point above z$n, the bin- wise technique effectively 
degenerates to a parameterized analysis at z > zsn- 



E. Nucleosynthesis 

In addition to the constraints on the expansion history 
(SNla, CMB, BAO and t ) we know that pde/ptot < 
0.045 (at 2cr confidence) during Big Bang Nucleosynthesis 
(BBN) [80] . Larger dark energy densities imply a higher 
expansion rate at that epoch (z ~ 6 x 10 8 ) which would 
result in a lower neutron to proton ratio, conflicting with 
the measured helium abundance, Y^ e . 



F. Dark Energy Parameterization 

Because of the variety of proposed dark energy mod- 
els, it has become usual to summarize observations by 
constraining a parameterized time-varying equation of 
state. Dark energy models are then confronted with ob- 
servations in this parameter space. The unique zeroth 
order parameterization of w is uu = wq (a constant), with 
w = — 1 characterizing the cosmological constant model. 
The observational data can be used to constrain the first 
derivative of w. This additional dimension in the DE 
parameter space may be useful in distinguishing models 
which have the same wq. From an observational stand- 
point, the obvious choice of 1st order parameterization 
is w(z) = w + ^z [8Xj |. This is rarely used today since 
currently considered DDE models are poorly portrayed 
by this functional form. The most popular parameteri- 



G. Summary of Current DDE Constraints 

If one assumes the popular wq — w a parameterization 
until last scattering, then all cosmological probes can be 
combined to constrain wq and w a . In a recent analysis 
of SNla, CMB and BAO observations, Davis et al. [681 ] 
found wo = — 1.0 ± 0.4 and w a = —0.4 ± 1.8 at 2a con- 
fidence (the contour is shown in Fig. [7]) . Using the same 
observations, Wood-Vasey et al. [591 ] assumed w a = and 
found w = w = -1.09 ± 0.16 (2a). 

The evolution of pde is related to w by covariant energy 
conservation [831 



= -3 (w(a) + 1) — . 

Pde a 



(7) 



The dark energy density corresponding to the wq — w a 
parameterization of w is thus given by 

Pde(z) = PdeO e 3 ^^ o -3(l+»0+™.). (g) 

The cosmic energy density history is illustrated in Fig. 
[3] Radiation and matter densities steadily decline as the 
dotted and dashed lines. With the DE equation of state 
parameterized as w(a) = wq + w a (l — a), its density his- 
tory is constrained to the light-grey area [68[. If the 
evolution of w is negligible, i.e. we condition on w„ k 0, 
then w(a) « Wq and the DE density history lies within 
the dark-grey band [59]. If the dark energy is pure vac- 
uum energy (or Einstein's cosmological constant) then 
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w = — 1 and its density history is given by the horizontal 
solid black line. 



IV. THE TEMPORAL DISTRIBUTION OF 
OBSERVERS 



in time. Thus our analysis will not depend on the nor- 
malization of this function and this uncertainty will not 
propagate into our analysis. There are also uncertainties 
in the location of the turnover at high redshift, and in 
the slope of the formation history at low redshift - both 
of these will affect our results. 



The energy densities p r , p m and pde , and the proximity 
parameter r we imagine we might have observed, depend 
on the distribution P t (t) from which we imagine our time 
of observation t b s has been drawn. What we can expect 
to observe must be restricted by the conditions necessary 
for our presence as observers [881 ] - Thus, for example, it 
is meaningless to suppose we might have lived during 
inflation, or during radiation domination, or before the 
first atoms 

We can, however, suppose that we are randomly se- 
lected cosmology-discovering observers, and we can ex- 
pect our observations of p m and pde to be typical of ob- 
servations made by such observers. This is Vilenkin's 
principle of mediocrity (9p| . Accordingly, the distribution 
Pt (t) for the time of observation t b s is proportional to the 
temporal distribution of cosmology-discovering observers 
(referred to henceforth as simply "observers"). Thus to 
solve the coincidence problem one must show that the 
proximity parameter we measure, r , is typical of those 
measured by other observers. 

The most abundant elements in the cosmos are hydro- 
gen, helium, oxygen and carbon [9l|. In the past decade 
> 200 extra solar planets have been observed via doppler, 
transit or microlensing methods. Extrapolation of cur- 
rent patterns in planet mass and orbital period are con- 
sistent with the idea that planetary systems like our own 
are common in the universe 92]. All this does not neces- 
sarily imply that observers are common, but it does sup- 
port the idea that terrestrial-planet-bound carbon-based 
observers, even if rare, may be the most common ob- 
servers. In the following estimation of Pt(t) we consider 
only observers bound to terrestrial planets. 
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FIG. 4: The terrestrial planet formation rate as estimated by 
Lineweaver 93] . It is based on a compilation of SFR measure- 
ments and has been corrected for the low metallicity of the 
early universe, which prevents the terrestrial planet forma- 
tion rate from rising as quickly as the stellar formation rate 
at z Z 4. 



The conversion from redshift to time depends on the 
particular cosmology, through the Friedmann equation, 



da\' 
dt) 



Hi 



H{afa 2 (10) 

Q r oa 2 + Q m o a 1 + 
n de0 exp[3w Q (a- 1)] a - 3w °- 3w '- 1 



In Fig.[5]we plot the PFR from Fig.[3]as a function of time 
assuming the best fit parameterized DDE cosmology. 



A. First the Planets... 

Lineweaver [§3| estimated the terrestrial planet forma- 
tion rate (PFR) by making a compilation of measure- 
ments of the cosmic star formation rate (SFR) and sup- 
pressing a fraction of the early stars f(t) to correct for 
the fact that the metallicity was too low for those early 
stars to host terrestrial planetary systems, 

PFR{t) = const x SFR(t) x f(t). (9) 

In Fig. |4] we plot the PFR reported by Lineweaver [93| 
as a function of redshift, z = - — 1. As illustrated in 

5 a 

the figure, there is large uncertainty in the normalization 
of the formation history. The number of stars orbited 
by terrestrial planets normalizes the distribution of ob- 
servers but, importantly, does not shift the distribution 



B. ... then First Observers 

After a star has formed, some non-trivial amount of 
time At i, s will pass before observers, if they arise at 
all, arise on an orbiting rocky planet. This time allows 
planets to form and cool and, possibly, biogenesis and the 
emergence observers. Ai b s is constrained to be shorter 
than the life of the host star. If we consider that our At b s 
has been drawn from a probability distribution PAt obs (t) ■ 
The observer formation rate (OFR) would then be given 
by the convolution 

/>oo 

OFR(t) = const x / PFR(T)P Atob ,(t - r)dr. (11) 
J o 

In practice we know very little about PAt ob3 (t) ■ It must 
be very nearly zero below about At b s ~ 0.5 Gyrs - this 
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FIG. 5: The terrestrial planet formation from Fig.[4]is shown 
here as a function of time. The transformation from redshift 
to time is cosmology dependent. To create this figure we 
have used best-fit values for the DDE parameters, wo — —1.0 
and w a = —0.4 68]. The y-axis is linear (c.f. the logarith- 
mic axis in Fig. [3} and the family of curves have been re- 
normalized to highlight the sources of uncertainty important 
for this analysis: uncertainty in the width of the function, 
and in the location of its peak. The observer formation rate 
(OFR) is calculated by shifting the planet formation rate by 
some amount At t s (= 4 Gyrs) to allow the planet to cool, 
and the possible emergence of observers. These distributions 
are closed by extrapolating exponentially in t. 



is the amount of time it takes for terrestrial planets to 
cool and the bombardment rate to slow down. Also, it 
is expected to be near zero above the lifetimes of sun- 
like stars (much above ~ 10 Gyrs). If we assume that 
our At bs is typical, then PAt ob3 (~t) has significant weight 
around At Q b s = 4 Gyrs - the amount of time it has taken 
for us to evolve here on Earth. 

A fiducial choice, where all observers emerge 4 Gyrs 
after the formation of their host planet, is PAt obs (t) = 
S(t— 4 Gyrs). This choice results in an OFR whose shape 
is the same as the PFR, but is shifted 4 Gyrs into the 
future, 

OFR(t) = const x PFR(t - 4 Gyrs) (12) 

(see the lower panel of Fig. [5]). Even for non-standard 
wq and w a values, this fiducial OFR aligns closely with 
the r(t) peak and the effect of a wider PAt obs is gen- 
erally to increase the severity of the coincidence prob- 
lem by spreading observers outside the r(t) peak. Hence 
using our fiducial PAt obs (which is the narrowest possi- 
bility) will lead to conclusions which are conservative in 
that they underestimate the severity of the cosmic coinci- 
dence. If another choice for PAt obs could be justified, the 



cosmic coincidence would be more severe than estimated 
here. We will discuss this choice in Section [VTl 

The OFR is then extrapolated into the future using a 
decaying exponential with respect to t (the dashed seg- 
ment in the lower panel of Fig. [5]). The observed SFH is 
consistent with a decaying exponential. We have tested 
other choices (linear & polynomial decay) and our results 
do not depend strongly on the shape of the extrapolating 
function used. 

The temporal distribution of observers P t (t) is propor- 
tional to the observer formation rate, 

P t (t) = const x OFR(t). (13) 

This observer distribution is similar to the one used by 
Garriga et al. [7^1 to treat the coincidence problem in a 
multiverse scenario. By comparison, our OFR(t) distri- 
bution starts later because we have considered the time 
required for the build up of metallicity, and because we 
have included an evolution stage of 4 Gyrs. Our distribu- 
tion also decays more quickly than theirs does. Some of 
our cosmologies suffer big-rip singularities in the future. 
In these cases we truncate Pt{t) at the big-rip. 

V. ANALYSIS AND RESULTS: DOES FITTING 
CONTEMPORARY CONSTRAINTS 
NECESSARILY SOLVE THE COSMIC 
COINCIDENCE? 

For a given model the proximity parameter observed 
by a typical observer is described by a probability distri- 
bution P r {r) calculated as 

p rM = E§W))- (14) 

The summation is over contributions from all solutions 
of t(r) (typically, any given value of r occurs at multiple 
times during the lifetime of the Universe). In Fig. [5] we 
plot P r {r) for the wq = —1.0, w a = —0.4 cosmology. In 
this case, observers are distributed over a wide range of 
r values, with 71% seeing r > r , and 29% seeing r < r . 

We define the severity S of the cosmic coincidence 
problem as the probability that a randomly selected ob- 
server measures a proximity parameter r lower than we 
do: 

S = P(r<r ) = l-P{r>r )= P r (r)dr. (15) 

Jr=Q 

For the wq = —1.0, w a = —0.4 cosmology of Figs. [5] 
and [13 the severity is S = 0.29 ± 0.09. This model does 
not suffer a coincidence problem since 29% of observers 
would sec r lower than we do. If the severity of the 
cosmic coincidence would be near 0.95 (0.997) in a par- 
ticular model, then that model would suffer a 2a (3a) 
coincidence problem and the value of r we observe really 
would be unexpectedly high. 

We calculated the severities S for cosmologies span- 
ning a large region of the wq — w a plane and show our 
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proximity r 



FIG. 6: The predicted distribution of observations of r is plot- 
ted for the parameterized DDE model which best-fits cosmo- 
logical observations: too = —1.0 and w a = —0.4. The prox- 
imity parameter we observe ro = ^ nSL ~ 0.35 is typical in 
this cosmology since only 29% of observers (vertical striped 
area) observe r < 0.35. The upper and lower limits on this 
value resulting from uncertainties in the SFR are 38% and 
20% respectively. Thus the severity of the cosmic coincidence 
in this model is S = 0.29 ± 0.09. This model does not suffer 
a coincidence problem. 



results in Fig. [7] using contours of equal S. The severity 
of the coincidence problem is low (e.g. S < 0.7) for most 
combinations of wq and w a shown. There is a coinci- 
dence problem, where the severity is high (S ?J 0.8), in 
two regions of this parameter space. These are indicated 
in Fig. [3 

Some features in Fig. are worth noting: 

• Dominating the left of the plot, the severity of the 
coincidence increases towards the bottom left-hand 
corner. This is because as wo and w a become more 
negative, the r peak becomes narrower, and is ob- 
served by fewer observers. 

• There is a strong vertical dipole of coincidence 
severity centered at (wq = 0, W a = 0). For (w 

0, w a > 0) there is a large coincidence problem be- 
cause in such models we would be currently wit- 
nessing the very closest approach between DE and 
matter, with pd e S> p m for all earlier and later times 
(see Fig. [5]:). For (wo ~ 0, w a < 0) there is an anti- 
coincidence problem because in those models we 
would be currently witnessing the DDE's furthest 
excursion from the matter density, with pd e and p m 
in closer proximity for all relevant earlier and later 
times, i.e., all times when Pt(t) is non-negligible. 

• There is a discontinuity in the contours running 
along w a = for phantom models (w < — 1). The 
distribution Pt(t) is truncated by big-rip singular- 
ities in strongly phantom models (provided they 
remain phantom; w a > 0). This truncation of late- 
time observers means that early observers who wit- 



ness large values of r represent a greater fraction of 
the total population. 

To illustrate these features, Fig. [5] shows the density 
histories and observer distributions for four specific ex- 
amples selected from the wq — w a plane of Fig. [7J 

We find that all observationally allowed combinations 
of wo and w a result in low severities (S < 0.4), i.e., there 
are large (> 60%) probabilities of observing the matter 
and vacuum density to be at least as close to each other 
as we observe them to be. 

It is not suggested that for arbitrary models P oos 
should be large when and only when pd e ~ p m - Indeed, it 
is easy to imagine pde ~ p m when there are not observers. 
Moreover, in some non-standard cosmological models it 
is possible to contrive pd e <C p m during times when ob- 
servers do exist. What our results suggest is that, for 
models parameterized with wq and w a satisfying current 
constraints, most observers (> 60%) will see pde and p m 
nearly equal (r > 0.35). 



VI. DISCUSSION 

Particle-theoretic approaches to the cosmic coincidence 
problem have focussed on the generation of a constant or 
slowly varying density ratio r. However it has not been 
made clear precisely how slowly the ratio r must evolve in 
order to solve the coincidence problem. In other words, 
the question "What DDE dynamics are required to solve 
the coincidence problem?" has not been addressed. 

In the present work we adopt the principle of medi- 
ocrity: that we should be typical observers, to try to 
answer this question. We estimate the temporal distri- 
bution of observers and devise a scheme for quantifying 
how unlikely the observation r > 0.35 is for an arbitrary 
DDE model. This scheme is applied to wo — w a paramet- 
ric models, and we identify regions of the wq — w a pa- 
rameter space in which the coincidence problem is most 
severe, however these are already strongly excluded by 
observations (see Fig. [TJ. 

Thus the main result of our analysis is that any realistic 
DDE model which can be parameterized as w — wq + 
w a (l — a) over a few e-folds, has pde ~ Pm for a significant 
fraction of observers. 

Central to our approach is the temporal distribution of 
observers as estimated using the distribution of terrestrial 
planets. Such observer selection effects are operating. 
Thus, while they may be difficult to quantify, they need 
to be considered whenever the cosmic coincidence is used 
to motivate new physics. These anthropic considerations 
operate in conjunction with (not in place of) fundamental 
explanations of the dark energy. 

Interacting quintessence models in which the proxim- 
ity parameter asymptotes to a constant at late times 

Ealilllllia.lailallMlaSliilSliilHl have 

been proposed as a solution to the coincidence problem. 
More recently, del Campo et al. (2^, [35[ have argued for 
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a broader class of interacting quintessence models that 
"soften" the coincidence problem by predicting a very 
slowly varying (though not constant) proximity parame- 
ter. Our analysis finds that r need not asymptote to a 
constant, nor evolve particularly slowly, partially under- 
mining the motivations for these interacting quintessence 
models. 

Caldwell et al. [6l[ and Scherrer [33[ have proposed 
that the coincidence problem may be solved by phan- 
tom models in which there is a future big-rip singularity 
because such cosmologies spend a significant fraction of 
their lifetimes in r ~ 1 states. In our work Pt(t) is termi- 
nated by big-rip singularities in ripping models. In non- 
ripping models, however, the distribution is effectively 
terminated by the declining star formation rate. There- 
fore the big-rip gives phantom models only a marginal 
advantage over other models. This marginal advantage 
manifests as the discontinuity along w a = on the left 
side of Fig. \T\ 

A running cosmological constant A(t) could arise from 
the renormalization group (RG) in quantum field the- 
ory [94], [9f| . The running lambda term can mimic 
quintessence or phantom behaviour and transit smoothly 
between the two |96[ . RG models represent interesting al- 
ternatives to scalar-field models of dark energy. In some 
variants [93] additional fields are introduced to address 
the cosmic coincidence problem by predicting a slowly 
varying density ratio r. Our results demotivate such ad- 
ditions and favor simplistic RG models. 

How strongly do these results depend on the assumed 
time it takes for observers to arise, At {, s ? In Lineweaver 
and Egan [121 ], where we performed an anlysis simi- 
lar to the present one (but limited to w — —1), we 
demonstrated that the results were robust to any choice 
A£ b s ~ [0, 11] Gyrs. However, for At b s > 12 Gyrs that 
analysis resulted in an unavoidable coincidence problem 
because most observers would arise late (during DE dom- 
ination) and would observe pde ^> Pm- The validity of the 
results of the present analysis are similarly limited. 

We could improve our analysis, in the sense of getting 
tighter coincidence constraints (larger severities), if we 
used a less conservative PAt obs - We used the most con- 
servative choice - a delta function - because the present 
understanding of the time it takes to evolve into observers 
is too poorly developed to motivate any other form of 
PAt ob3 ■ Another possible improvement is the DE equa- 
tion of state parameterization. We used the current stan- 
dard, w — wo + w a (l — a), which may not parameterize 
some models well for very small or very large values of a. 

We conclude that DDE models need not be fitted 
with exact tracking or oscillatory behaviors specifically to 
solve the coincidence by generating long or repeated pe- 
riods of pde ~ Pm- Also, particular interactions guaran- 
teeing pde ~ Pm for long periods are not well motivated. 
Moreover phantom models have no significant advantage 
over other DDE models with respect to the coincidence 
problem discussed here. 
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FIG. 7: Here we plot contours of equal severity S in wo — w a parameter space. S is the fraction of observers who see r < ro- If 
S is large, a large percentage of observers should see r lower than we do - those models suffer coincidence problems. The thick 
black contour represents the observational constraints on wo and w a from Davis et al. [68| ] (2a confidence and marginalized over 
other uncertainties). In Lineweaver and Egan pl| we showed that the severity of the coincidence problem is low for ACDM 
(indicated by the "+"). Values of wo and w a that result in a mild coincidence problem (e.g. S i5 0.7) are already strongly 
excluded by observations. This leads to our main result: none of the models in the observationally allowed regime suffer a 
cosmic coincidence problem when our estimate of the temporal distribution of observers P b s (t) is used as a selection function. 
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scalefactor a 

FIG. 8: History of the energy densities in radiation (dotted line), matter (dashed line) and dark energy (thick black line) for 
four parameterized DE models from Fig. [7] The proximity parameter r (thin black line) and the temporal distribution of 
observers Pt(t) (grey shade) are also given. Panel a shows a phantom model with a constant equation of state w = —3.5. In 
this model the phantom density increases quickly and the r(t) peak is narrow. As a result, a large fraction of observers live 
while the matter and dark energy densities are vastly different (r « 0) and there is a mild coincidence problem (S « 0.8). 
This might be used to rule-out the model shown in Panel a, except that it is already strongly excluded by direct cosmological 
observations (refer to Fig. [TJl . Panel b shows a phantom model which lies within the observationally allowed 2er region. There 
is no coincidence problem in this model (S w 0.4). Panel c shows a model in which there is a coincidence problem (S w 0.95). 
This models lies within the cluster of contours in the upper right-hand corner of Fig. [7] In this model the dark energy dominates 
the past and future energy budget. Again however, the coincidence problem can tell us nothing new, as this model is already 
strongly excluded by observations. Panel d shows a model in which there is an anti-coincidence problem. This models lies 
within the cluster of contours in the lower right-hand corner of Fig. [7] In this model the dark energy and matter densities are 
more similar (r is greater) in the recent past and near future (although r — > further into the past or future). According to 
the observer distribution Pt(t) most observers live near the current epoch, during r > 0.35, with just 7% living during r < 0.35 
(S — 0.07) in this particular model. One might argue that this model can be ruled out because our value of r is anomalously 
small. However, this model too is already strongly excluded by observations. 
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APPENDIX A: NUMERICAL VALUES FOR 
PARAMETERS OF MODELS ILLUSTRATED IN 
FIG. [2] 



TABLE I: Free parameters of the DDE models illustrated in 
Fig. [2] These values were chosen such that observational con- 
straints are crudely satisfied. These are by no means the only 
combinations fitting observations. These values are intended 
for the purposes of illustration in Fig. [2] Units are Planck 
units. 

Model Parameter Value 



power law tracker quintessence 


a 


2 








M 


1.4 x 


10" 


124 


exponential tracker quintessence M 


1.3 x 


10" 


124 


tracking oscillating energy 


M 


1.8 x 


10" 


126 




X 


4 








A 


0.99 








V 


2.7 






interacting quintessence 


A 


1.4 x 


10" 


119 




B 


9.7 








C 


10 






Chaplygin gas 


a 


1 








A 


2.8 x 


10" 


l>4<; 
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